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Abstract. This paper develops the notion of implicit Lagrangian systems on 
Lie algebroids and a Hamilton-Jacobi theory for this type of system. The Lie 
algebroid framework provides a natural generalization of classical tangent bun- 
dle geometry. We define the notion of an implicit Lagrangian system on a Lie 
algebroid E using Dirac structures on the Lie algebroid prolongation 7 E E* . 
This setting includes degenerate Lagrangian systems with nonholonomic con- 
straints on Lie algebroids. 
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1. Introduction 

There is a vast literature on the Lagrangian formalism in mechanics, which 
is due to the central role played by these systems in the foundations of modern 
mathematics and physics. In many interesting systems, problems often arise due 
to their singular nature, which gives rise to constraints that address the fact that 
the evolution problem is not well-posed (internal constraints) . Constraints can also 
manifest a priori restrictions on the states of the system which arise due to physical 
arguments or from external conditions (external constraints). Both cases are of 
considerable importance. 

Systems with internal constraints are quite interesting since many dynamical 
systems are given in terms of presymplectic forms instead of the more habitual 
symplectic ones. The more frequent case appears in the Lagrangian formalism of 
singular mechanical systems which are commonplace in many physical theories (as 
in Yang-Mills theories, gravitation, etc). 

Systems subjected to external constraints (holonomic and nonholonomic) have 
a wide range of applications in many different areas: engineering, optimal control 
theory, mathematical economics (growth economic theory), subriemannian geom- 
etry, motion of microorganisms, etc. Interconnected and implicit systems play a 
key role in, for example, controlled mechanical systems like robots. An impor- 
tant class of implicit mechanical systems is those with nonholonomic constraints, 
which has a long and rich history (see, for instance, [4] and [35]). The Lagrangian 
and Hamiltonian approaches for such systems have been extensively developed (see 
[13 [MIES EI]), including symmetry and reduction (see [51 15] |22] [M] 1213] ) . 

Some authors have given descriptions of L-C circuits and nonholonomics systems 
in the context of Poisson structures (see [55J[5B]) an d later in the general context of 
Dirac structures (see [5] [35]) from a Hamiltonian point of view. Inspired by these 
works, Yoshimura and Marsden in [JSJ HI] have developed a Lagrangian formalism 
making use of the framework of Dirac structures. 

Recent investigations have lead to a unifying geometric framework covering a 
plethora of particular situations. It is precisely the underlying structure of a Lie 
algebroid on the phase space which allows a unified treatment. This idea was first 
introduced by Weinstein |42) in order to define a Lagrangian formalism which is 
general enough to account for different types of systems. The geometry and dy- 
namics on Lie algebroids have been extensively studied during the past years. In 
particular, in [30], E. Martinez developed a geometric formalism of mechanics on 
Lie algebroids similar to Klein's formalism of ordinary Lagrangian mechanics and, 
more recently, a description of the Hamiltonian dynamics on a Lie algebroid was 
given in [271 [3T]. The key concept in this theory is the prolongation, r J E E, of the 
Lie algebroid over the fiber projection r (for the Lagrangian formalism) and the 
prolongation, 7 E E*, over the dual fiber projection r* : E* — > Q (for the Hamil- 
tonian formalism). See [27 for more details. Of course, when the Lie algebroid 
is E = TQ we obtain that 7 E E = T(TQ) and 7 E E* = T(T*Q), recovering the 
classical case. Another approach to the theory was discussed in [T7] . 

The notion of nonholonomic systems on a Lie algebroid was introduced in [9] 
when studying mechanical control systems and an approach to mechanical systems 
on Lie algebroids subject to linear constraints was presented in [34] , A recent 
comprehensive treatment of nonholonomic systems on a Lie algebroid has been 
develop in [10| , where the authors identify suitable conditions guaranteeing that the 
system admits a unique solution and show that many of the properties that standard 
nonholonomic systems enjoy have counterparts in the Lie algebroid setting. 
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On the other hand, singular or degenerate Lagrangian systems and vakonomic 
mechanics on Lie algebroids (obtained through the application of a constrained vari- 
ational principle) also have been studied. In |21j . the authors introduce a constraint 
algorithm for presymplectic Lie algebroids which generalizes the well-known Gotay- 
Nester- Hinds algorithm (see jTS] ) and applies it to singular Lagrangian systems on 
Lie algebroids. Moreover, they develop a geometric description of vakonomic me- 
chanics on Lie algebroids using again the constraint algorithm. 

As a consequence of all these investigations, one deduces that there are several 
reasons for discussing unconstrained (constrained) Mechanics on Lie algebroids: 

i) The inclusive nature of the Lie algebroid framework. In fact, under the same 
umbrella, one can consider standard unconstrained (constrained) mechanical sys- 
tems, (nonholonomic and vakonomic) Lagrangian systems on Lie algebras, uncon- 
strained (constrained) systems evolving on semidirect products or (nonholonomic 
and vakonomic) Lagrangian systems with symmetries. 

ii) The reduction of a (nonholonomic or vakonomic) mechanical system on a 
Lie algebroid is a (nonholonomic or vakonomic) mechanical system on a Lie alge- 
broid. However, the reduction of a standard unconstrained (constrained) system 
on the tangent (cotangent) bundle of the configuration manifold is not, in general, 
a standard unconstrained (constrained) system. 

iii) The theory of Lie algebroids gives a natural interpretation of the use of quasi- 
coordinates (velocities) in Mechanics (particularly, in nonholonomic and vakonomic 
mechanics). 

On the other hand, Hamilton-Jacobi theory has been studied for different type of 
systems for many years. For degenerate Lagrangian systems, some work have been 
done on extending Hamilton-Jacobi theory, using Dirac's theory of constraints (see, 
e.g., [18]) and from a geometric point of view (see [7]). For nonholonomic systems, 
in [20], Iglesias-Ponte, de Leon and Martin de Diego generalized the geometric 
Hamilton-Jacobi theorem (see Theorem 5.2.4. in to nonholonomic systems, 
which has been studied further (see ;8, 36, 37 ). More recently, in [25], the authors 
have presented a Hamilton-Jacobi theory which can deal with both degeneracy and 
nonholonomic constraints. In the context of Lie algebroids, de Leon, Marrero and 
Martin de Diego have developed a more general formalism which is also valid for 
for nonholonomic systems on a Lie algebroid (see [H]), and, in [2J, the authors 
have presented a Hamilton-Jacobi equation for a Hamiltonian system on a skew- 
symmetric algebroid. 

The goal of this paper is to generalize Hamilton-Jacobi theory to implicit La- 
grangian systems on a Lie algebroid based on Dirac structures. We introduce the 
notion of an implicit Lagrangian system on a Lie algebroid E using the induced 
generalized Dirac structure D\i on the Lie algebroid prolongation 7 E E* that is nat- 
urally induced by a vector subbundle 11 of E and we obtain the Hamilton-Jacobi 
theorem for this kind of systems. This setting includes degenerate Lagrangian 
systems with nonholonomic constraints. 

The paper is organized as follows. In Section [2] we collect some preliminary 
notions and geometric objects on Lie algebroids, including differential calculus, 
morphism and prolongations. We also recall the definition and some properties of 
(generalized) Dirac structures on vector spaces, vector bundles and manifolds. In 
Section|3l first we introduce and study the generalized Dirac structure T>u on 7 E E* 
induced by a vector subbundle U of the Lie algebroid E. The main goal of this 
section is to define implicit Lagrangian systems in terms of induced Dirac structures. 
In Section U] we develop a Hamilton-Jacobi theory for implicit Lagrangian systems 
on a Lie algebroid. We apply the results obtained to some particular cases, in 
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Section [SJ recovering some known results. The paper ends with our conclusions 
and a description of future research directions. 

2. Preliminaries 

2.1. Lie algebroids. Let £ be a vector bundle of rank n over a manifold Q of 
dimension m and r : E — > Q be the vector bundle projection. Denote by T(E) the 
C°°((5)-module of sections of r : i? — >• Q. A iie algebroid structure ([•, •], p) on 
i? is a Lie bracket [•, •] on the space and a bundle map p : E — > TQ, called the 
anchor map, such that if we also denote by p : Y{E) — >■ X(Q) the homomorphism 
of C°°((5)-modules induced by the anchor map, then 

lXJY] = flX,Y}+p(X)(f)Y, 

for X,Y 6 T(E) and / G C°°(Q). The triple (£, [•, -],p) is called a Lie algebroid 
over Q (see [2S])- 

If (E, [•, -],/o) is a Lie algebroid over Q, then the anchor map p : T(E) — > X(Q) 
is a homomorphism between the Lie algebras (T(E), [•, •]) and (X(Q), [■, ■]). 

Standard examples of Lie algebroids are real Lie algebras of finite dimension and 
the tangent bundle TQ of an arbitrary manifold Q. In more detail, let (g, [■, -] ) 
be a real Lie algebra of finite dimension. Then, consider the vector bundle r : 
g — > { one point }. The section of this vector bundle can be identified with the 
elements of g and, therefore, we can consider the Lie bracket given by the Lie 
algebra structure [-,-] B on q and the anchor map p given by the null map. So, 
(fl) [•) 'IbjO) is a Lie algebroid over a point. On the other hand, let Q a manifold. 
The sections of the tangent bundle te — tq : E = TQ — >■ Q may be identified with 
the vector fields on Q, the Lie bracket on T(te) = X(Q) is the usual vector fields 
bracket and the anchor map is the identity on TQ. Then, the triple (TQ, [•, ■], Id) 
is a Lie algebroid over Q. 

Another example of a Lie algebroid may be constructed as follows. Let ir : P — > Q 
be a principal bundle with structure group G. Denote by (f> : G x P — > P the free 
action of G on P and by TQ : G x TP TP the tangent lifted action of G on TP. 
Then, one may consider the quotient vector bundle Tp\G : TP/G — > Q = P/G and 
the sections of this vector bundle may be identified with the vector fields on P which 
are invariant under the action $. Using the fact that every G- invariant vector field 
on P is 7r-projectable and the fact that the standard Lie bracket on vector fields 
is closed with respect to G- invariant vector fields, we can induce a Lie algebroid 
structure on TP/G. The resultant Lie algebroid is called the Atiyah (gauge) 
algebroid associated with the principal bundle ir : P — > Q (see [2"T1 128j ). 

Now, let (E, [•,•], p) be a Lie algebroid, then one may define the differential 

of E, d E : T{A k E*) -> T(A fe+1 £;*), as follows 

fc 

d E p(X , ...,X k ) = ]T(-l)VPQ)(/i(A , V X k )) 

i=0 

+ Xj], x , v V,,... .,X k ), 

i<j 

for u G r(A fc £:*) and X Q ,...,X k e T(E). It follows that (d E ) 2 = 0. Moreover, if 
X G r(i?), one may introduce, in a natural way, the Lie derivative with respect 
to X, as the operator £ E : T(A k E*) -)■ T{A k E*) given by £f- = i x od E + d E oi x . 

Note that if E = TQ and X G T(E) = X(Q) then d T ® and £ X Q are the usual 
differential and the usual Lie derivative with respect to X, respectively. 



If we take local coordinates (x l ) on an open subset U of Q and a local basis {e Q } 
of sections of E defined on U, then we have the corresponding local coordinates 
(x t ,y a ) on E, where y a (e) is the a-th coordinate of e G £7 in the given basis. Such 
coordinates determine local functions p l a , C 7 ^ on Q which contain local information 
about the Lie algebroid structure and, accordingly, they are called the structure 
functions of the Lie algebroid. They are given by 



D7 e 



and p(e c 



Pn 



dx l 



These functions should satisfy the relations 

d Pg 



pi 



cyclic(ot,ft ,7) 



dx l 



— P~f^a/3> 



= 0, 



which are usually called i/ie structure equations. 

If / € C°°(Q), we have that 

where {e a } is the dual basis of {e a }. On the other hand, if 6 G T(E*) and 6* = 6* 7 e 7 , 
it follows that 



(2.1) 



Ae 7 . 



In particular, 



■-C^Ae 7 . 



2.2. Morphisms. Let [•, •], p) and (£", [•, •]', //) be Lie algebroids over Q and 
Q', respectively. A morphism of vector bundles (F, /) from E to E' 



E 



Q 



f 



E' 



Q' 



is a Lie algebroid morphism if 

d E {{FJ)*<t> , ) = {F,fy{d E ' ( l>>), for 0'GL(A fc (£')*). (2.2) 
Note that (F, f )*</>' is the section of the vector bundle A k E* — > Q defined by 

({F, /)V)«(<*i, ■ • • , Ofc) = 0/ ( ,)(^(ai), ■ • ■ , ^M), 

for x € Q and 01, . . . , G E x , where E x denotes the fiber of E at the point 16Q. 
We remark that l|2.2[) holds if and only if 

° /) = {F, fT(d E 'g'), for G C°°(Q')> 
/)*«') = (F,/)*(d £ 'a'), for a' G r((£')*)- 

If (F, f) is a Lie algebroid morphism, / is an injective immersion and F\E m : E x — > 
Ef(x) ^ s injective, for all x G Q, then (U, [•, •], p) is said to be a 2/ie subalgebroid 

of (£',[•, •]>')■ 



6 



M. LEOK AND D. SOSA 



If Q = Q' and / = id : Q — > Q then, it is easy prove that the pair (F, id) is a Lie 
algebroid morphism if and only if 

F{X, Y] = {FX, FY]', p'(FX) = p(X), 

for X,Y e T(E). 

2.3. Poisson structure on E* . Let (E, [•, -},p) be a Lie algebroid over Q and E* 
be the dual bundle to E. Then, E* admits a linear Poisson structure that is, 
He* is a 2-vector on E* such that 

[IIe*,IIe*] = 0, 

and if / and /' are linear functions on E* , we have that He* {d TE /, d TE /') is also 
a linear function on E* . If (x l ) are local coordinates on Q, {e a } is a local basis 
of T(E) and (x l ,p a ) are the corresponding local coordinates on E* , then the local 
expression for LT^* is 

where p l a and are the structure functions of E with respect to the coordinates 
(V) and to the basis {e Q }. The Poisson structure ILs* induces a linear Poisson 
bracket of functions on E* which we will denote by { , }e*- In fact, if F, G £ 
C°°(E*) then 

{F,G} E * =U E *(d TE 'F,d TE 'G). (2.3) 
(For more details, see [2"T]). 

2.4. The prolongation of a Lie algebroid over a flbration. Let (E, [•, •], p) 

be a Lie algebroid of rank n over a manifold Q of dimension m and tt : P — > Q be 
a flbration, that is, a surjective submersion. 

We consider the subset J E P of Ex TP defined by J E P = [J 7 E P, where 

veP 

7 E P = {(&, w) e ^ (p) x T p P | p(b) = (T p tt)(v)}, 
and T7r : TP — > TQ is the tangent map to 7r. 
Denote by r 71 " : T B P — )■ P the map given by 

T«(b,v)=T P {v), 

for (b,v) G T B P, where Tp : TP — > P is the canonical projection. Then, if m! is 
the dimension of P, one may prove that 

dim 7 E P = n + w' — m. 

Thus, we conclude that T B P is a vector bundle over P of rank n + m! — m with 
the vector bundle projection t" : 7 E P — > P. 

A section X of t 1 : T^P — >• P is said to be projectable if there exists a section 
X of r : B — > Q and a vector field U on P which is 7r-projectable to the vector field 
and such that X(p) = (X(ir(p)), U(p)), for all p G P. For such a projectable 
section X, we will use the following notation X = (X, U). It is easy to prove that 
one may choose a local basis of projectable sections of the space T(T B P). 

The vector bundle t w : 7 E P — > P admits a Lie algebroid structure ([•, •] 7r , p 7r ). 
In fact, 

[(XuUt), (X 2 , u 2 )Y - (lX 1 ,X 2 j, [U u U 2 ]), p'iXuU!) = Ul 

The Lie algebroid (7 E P, [•, -J^, p*) is called the prolongation of E over tt or the 
E -tangent bundle to P. Note that if pr x : 7 E P — » E is the canonical projection 
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on the first factor, then the pair (pr 1; 7r) is a morphism between the Lie algebroids 
(7 E P, [•, ■j 7r ,p 7T ) and (E, [•, -],p) (for more details, see [37]). 

Example 2.1. Let (E, [•, be a Lie algebroid of rank n over a manifold Q 
of dimension m and r : E — >• Q be the vector bundle projection. Consider the 
prolongation 7 E E of 2£ over r, 

T £ £7 = {(e,w) e E x TE\ p(e) = (Tr)(u)}. 

T B _E is a Lie algebroid over E of rank 2n with Lie algebroid structure ([•, -] r , p T ). 

If (x 1 ) are local coordinates on an open subset U of Q and {e a } is a basis of 
sections of the vector bundle r _1 ([/) — > U, then {X Q , V a } is a basis of sections of 
the vector bundle (r T )" 1 (r" 1 (C/)) ->■ r _1 (£/"), where t t : T B P ->■ £ is the vector 
bundle projection and 

I a (a) = (e a (r(e)),p>Aj, V a (e) = (o, ^ J, (2.4) 

for e £ t _1 ([/). Here, are the components of the anchor map with respect to 
the basis {eo,} and (x l ,y a ) are the local coordinates on E induced by the local 
coordinates (x l ) and the basis {e a }. Using the local basis {XcVq,}, one may 
introduce, in a natural way, local coordinates (x l ,y a ; s a , w a ) on 7 E E. If cu is a point 
of (r T ) _1 (r _1 ([/)), then (x l ,y a ) are the coordinates of the point t t {uj) e r _1 (J7) 
and 

u = s a X a {T T (uj)) + w a V a (T T (uj)). 
On the other hand, we have that 

[x„, x p y = e^x 7 , [x Q , v P Y = lv a ,v ] T = o, 

for all a and /3, where are the structure functions of the Lie bracket [•, •] with 
respect to the basis {e Q }. 

The vector subbundle (7 E E) V of 7 E E whose fiber at the point e e E is 
{7 E E) V = {(0,v) eEx T e E\(T e r)(v) = 0} 

is called the vertical subbundle. Note that (7 E E) V is locally generated by the 
sections {V Q }. 

Two canonical objects on 7 E E are the Euler section A and the vertical 
endomorphism S. A is the section of 7 E E — > E locally defined by 

A = y a V a , 

and S is the section of the vector bundle (7 E E) (£> (7 E E)* — »■ E locally characterized 
by the following conditions 

S(X a )=V a , S(V a )=0, for all a. (2.5) 

Finally, a section £ of 7 E E — > E is said to be a second-order differential equa- 
tion (SODE) on E if = A or, alternatively, pr x (£(e)) = e, for all e e E (for 
more details, see [37J ). 

Example 2.2. Let (I?, [•, •], p) be a Lie algebroid of rank n over a manifold Q of 
dimension m and r* : £?* — >■ Q be the vector bundle projection of the dual bundle 
E* to E. 

We consider the prolongation 7 E E* of £7 over t* , 

7 E E* = {(e', v) e E x TE* \ p(e') = (Tt*){v)}. 
7 E E* is a Lie algebroid over E* of rank 2n with Lie algebroid structure ([•, -] r , p T ). 
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If (x l ) are local coordinates on an open subset U of Q, {e a } is a basis of sections 
of the vector bundle r _1 (£/) — > U and {e a } is the dual basis of {e Q }, then {^ a , T a } 
is a basis of sections of the vector bundle (t t )~ 1 ((t*)~ 1 (U)) — > (r*) _1 (C/), where 
t t : 7 E E* — > E 1 * is the vector bundle projection and 

MO-^O^J, ^( e *) = (o,Aj, (2 . 6) 

for e* 6 (r*) _1 (i7). Here, (x l ,p Q ) are the local coordinates on E* induced by the 
local coordinates (x l ) and the basis {e Q } of T(E*). Using the local basis {^ a , V*}, 
one may introduce, in a natural way, local coordinates (x l ,p a ; z a , u a ) on 7 E E* . If 
uj* is a point of (r T )~ 1 ((t*)~ 1 (U)), then (x l ,p Q ) are the coordinates of the point 
t t *(w*) e (r*)- x (l7) and 

^* = 2 Q y Q (r r * («*)) + UQ y Q (r T * (w*)). 
On the other hand, we have that 

IMa, ^] T * = e^y 7 , [y Q) ?T * = pp a , &Y* = o, 

■ . 8 (2-7) 

for all a and /3. Thus, if {T^a} is the dual basis of {V then 



3> 7 = 0, 



for / G C°°(E*). 

We may introduce a canonical section X E of the vector bundle (7 E E*)* — > _E* 
as follows. If e* e E* and (e, u) is a point of the fiber of 7 E E* over e*, then 

X E (e*)(e,v) = (e*,e), 

where (•, •) denotes the natural pairing between E* and £\ As is called the Liou- 
ville section of (T E E*)*. 

Now, i/ie canonical symplectic section £l E is the nondegenerate closed 2- 
section defined by 

Then, we have that the map Q E : T E E* -> (T B ^*)* defined as 

& E (X)=i x Q E , (2.8) 

for all X G 7 E E*, where ix denote the contraction by X, is a vector bundles 
isomorphism. 

In local coordinates, 

\ E {x\p a ) = p a y a , 

tt E (x\ Pa ) = TA? a + \& aP Pji a A (2.9) 

Remark 2.3. The linear Poisson bracket {-,-} E * on E* induced by the Lie al- 
gebroid structure on E (see l|2.3p ) can be also defined in terms of the canonical 
symplectic 2-section fl E . In fact, for F,G G C°°(E*), we have that 

{F,G} E , = n E ((n E )-\d 7BE *F),(rt E )-\d' xEE *G)). 
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2.5. Dirac structures. In this section we briefly recall the definition and some 
properties of Dirac structures on vector spaces, vector bundles and manifolds (see 
[TT1 IT2"]). The construction of a Dirac structure will be reviewed, which will be 
important for defining implicit Lagrangian systems. 

Let V be an n-dimensional vector space, V* be its dual space, and let (•, •) be 
the natural pairing between V* and V. A Dirac structure on V is a subspace 
D C V © V* such that D = D ± , where D 1 - is the orthogonal complement of D, 
that is, 

D x = {(«, (3)eV®V*\ (/3, v) + (a, u) = 0, for all (v, a) G D}. 

It is easy to prove that a vector subspace D C y © is a Dirac structure on V 
if and only if dim D = n and (a, v) + (a, u) = 0, for all (v, a), (v, a) € 2). From the 
definition of a Dirac structure, for each (v, a) G 2), we have that (a, u) = 0. 

If V is a vector bundle over a manifold Q, let V ©q y* be the Whitney sum 
bundle over Q, that is, it is the bundle over the base Q and with fiber over the point 
x G Q equal to V x x V*, where V x (respectively, V*) is the fiber of V (respectively, 
V*) at the point x. A Dirac structure on V is a subbundle D C V (Bq V* that 
is a Dirac structure in the sense of vector spaces at each point x G Q. 

Now, let M be a smooth diffcrentiable manifold and tm : TM — > M its tangent 
bundle. An almost (in the terminology of [53]) or generalized (in the terminology 
of [13 ) Dirac structure on M is a subbundle D C TM @m T*M which is a Dirac 
structure in the sense of vector bundles. 

In geometric mechanics, almost Dirac structures provide a simultaneous gener- 
alization of both 2-forms (not necessarily closed and possibly degenerate) as well 
as almost Poisson structures (that is, bracket that need not satisfy the Jacobi iden- 
tity). A Dirac structure on M is an almost Dirac structure that additionally 
satisfies the following integrability condition 

(£ Xl a2,X z ) + (£x 2 a 3 ,Xi) + {£x 3 ai,X 2 ) = 0, 

for all (Xi, ai), (X 2 , a^), (X 3 , 013) £ D, and where £x denotes the usual Lie de- 
rivative with respect to the vector field X. This generalizes closedness for the 
symplectic form, and the Jacobi identity for Poisson structures. For the remainder 
of this paper, we will primarily be concerned with almost Dirac structures, since it 
allows one to incorporate nonholonomic constraints. 

Two constructions of almost Dirac structures on a manifold are given as follows. 
The first construction is induced by a distribution and a 2-form on the manifold. 
Let M be a manifold, £1 be a 2-form on M and Am be a distribution on M . Denote 
by £! b the associated flat map and by C T*M the annihilator of Am- Then, 
from Theorem 2.3 in [43], we have that D M C TM © M T*M defined, for each 
x e M, by 

?>m{x) ={{v x ,a x ) £ T X M x T*M \ v x G A M (x) and 
a w -n\x)(v x )eA° M (x)} 

is an almost Dirac structure on M (see also Theorem 3.2 in [T5]). 

For the case when M = T*Q and 57 = Ht*q is the canonical symplectic 2-form, 
this almost Dirac structure was used to introduce the notion of implicit Lagrangian 
systems in standard mechanics (see [4*51 FH] ). 

The second almost Dirac structure is induced by a codistribution and a skew- 
symmetric 2-tensor on the manifold. Let M a manifold, IT : T*M x T*M ->Rbea 
skew-symmetric 2-tensor and A M a codistribution on M . Denote by j)n : T*M — > 
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TM the associated sharp map and by ker A* M the distribution on M defined as 

ker A* M = {X £ TM \ a(X) = 0, for all a £ A* M }. 

Then, from Theorem 2.4 in [13], we have that D M C TM © M T*M defined, for 
each x £ M, by 

BmW ={(«». Ox) € T a M x T x *Af I a x £ A* M {x) and 

- $n(x)(a x ) £ ker A^(a;)} 

is an almost Dirac structure on M (see also Theorem 3.1 in [T3"]). 

For the case when M = T*Q and II = IIt*q is the canonical Poisson structure 
on T*Q, this almost Dirac structure coincides with the almost Dirac structure 
described before which was used to introduce the notion of implicit Lagrangian 
systems in standard mechanics (see [4"31 HI] ). 

3. Implicit Lagrangian systems on a Lie algebroid 

3.1. Induced almost Dirac structure. First, we introduce the notion of an 
induced almost Dirac structure on the Lie algebroid prolongation 7 E E* of a 
Lie algebroid r : E — > Q. This almost Dirac structure is induced by a vector 
subbundle XL of E, that is, U C E such that tu — rm : U — > Q is a vector bundle. 

Consider the dual vector bundle r* : E* — > Q of r : E — > Q. We can define its 
prolongation to the corresponding prolongation Lie algebroids Tr* : 7 E E* — > 7 E Q 
as the identity in the first component and the tangent map of r* in the second, that 
is, Tr* = (id, Tr*). It is easy to prove that it is a Lie algebroid morphism between 
7 E E* — > E* and 7 E Q — > Q (see [32] for a general definition of the prolongation of a 
map). Moreover, we can identify 7 E Q with E and then Tr* = pr x , pr x : 7 E E* — > E 
being the projection on the first factor. 

The vector subbundle U can be lifted to a vector subbundle U<je e , C 7 e E* as 
follows 

U 7EE , = (p^)- 1 ^). (3.1) 
Denote by U° B£ ,, C (7 E E*)* its annihilator. Then, we have the following result. 

Theorem 3.1. Let (E, [•, ■],/?) be a Lie algebroid over a manifold Q and 11 be a 
vector subbundle of E. For each e* £ E* , let 

V u (e*) ={(X e »,a e *) £ 7 E E* x (7 E E*)*\X e , £ U 7 B E ,(e*) and 

a e * -n E (e*)(X e .) €U° EE ,(e*)}. 

Then, T> u C 7 E E* ® E , (7 E E*)* is an almost Dirac structure on 7 E E* . 

Proof. First, it is not difficult to prove that, since \L<je e * is a vector subbundle of 
7 E E*, T> u is a vector subbundle of 7 E E* © B . (7 E E*)*. 

Second, the orthogonal of D u C 7 E E* @e* (7 e E*)* is given at e* £ E* by 
K(e*) ={(Y e ,,(3 e ,) £ 7 E E* x (7 E E*)* \ a e *(Y e *) + /3 e *(X e .) - 0, 

for all X e . £ U 7 E E ,{e*) and a e * - VL E {e*){X e .) £ U° 7 EE . (e*)}. 

To check that D u (e*) C V^(e*), we consider (X e *,a e *) £ D u (e*) and then, for 
any (X' e ,,a' e ,) £ D u (e*), we have that 

a e .{X' e ,) + a' e «(X e „) = n E (e*)(X e .,X' e «) + n E (e*)(X' e ,,X e .) = 0, 

by the skew-symmetry of £l E . This implies that (X e »,a e ») £ D^(e*). Therefore, 

Vu(e*) C Di(e*). (3.3) 
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Now, to prove that T>£(e*) C D u (e*), let (F e «,/3 e ») G £>u( e *)- Then, we have 
that 

a e *{Y e .) + p e ,{X e ») = 0, (3.4) 
for all (X e »,a e ») G Tf,£* x {7 E E*)* such that X e . G Ho-^ (e*) and a e « - 
n^(e*)(X e .) G U° 7EE ,(e*). If we choose X e * = and a e . G lX^ ES »(e*), then 
(X e *,Q! e *) G I>u(e*). Therefore, using (|3.4p . we obtain that a e *(Y" e *) = 0, for all 
a e * G U^ B£ ,„(e*). Then, we conclude Y e * G U T B£;* (e*). On the other hand, let 
X e * G \i'jE E »{e*) be arbitrary and suppose that a e »(Z e .) = ^£;(e*)(X e » , Z e . ), for 
allZ e » eU 7 E E ,(e*). Since y e . G U 7 e e , (e*), we have a e * (Y e . ) = Q E (e*)(X e , , Y e . ) 
and, from (|3 .4[) . we deduce that 

il E (e*)(X e *,Y e *) + p e *(X e *) = 0, 

for all X e . G U T E E ,(e*). This implies that /3 e . -n^(Y e ») G U° B£ , (e*). Therefore, 
(Y e *,/3 e *) G Du(e*) and thus 

2) l7 (e*) c D u (e*). (3.5) 

Given (|3.3[) and (|3.5j) . we conclude that D^(e*) = T> u (e*), and the result follows. 

□ 

In what follows, we will obtain a local representation of the almost Dirac struc- 
ture T>u induced on 7 E E* by a vector subbundle U of E. Consider local coordinates 
(x % ) on Q, a local basis {e a } of sections of E and the corresponding local coordi- 
nates {x\y a ) on E. Let {^ a , IP"} be the local basis of r r * : ->■ defined 
by (|2.6[) induced by the local coordinates on Q and the local basis {e a } of E 1 
and (x' t ,p a ; z a ,u a ) be the induced local coordinates on 7 E E* . 

Thus, we can locally represent the fiber of H<je e * at a point (x l ,p a ) G E* as 

U M (x> Q ) = {(x\p a ;z a ,u a )\(x\z a ) 6 11(4 

If we denote by (a;*,^; r Q , u Q ) the corresponding local coordinates induced on 
(T E E*)* by the dual basis {y a ,9 a } of 3> Q }, then the annihilator of U t e e * 
is locally given by 

U° 7 E E ,{x\p a ) = {(x\p a ;r a ,v a )\v a = and (x\r a ) G U°(x 1 )}. 
From (|2.9[) . we have that 

f4(:c\p Q )(x\p Q ;z Q ,u Q ) = (a;',p Q ; -u a - Q^^z 13 ,z a ) (3.6) 
and then the condition a e * — rt b E (e*)(X e * ) G U^. Ef; »(e*) can be written locally as 

v a =z a and {x\r a +u a + Ql p p^z p ) eU°(i'), 
where X e * = (x\p a ; z a , u a ) and a e » = (ir^Pa,; r Q , u Q ). 
Finally, we obtain that 

T> u (e*) ={{X e *,a e .) G 7 E E* x (7 E E*)* | X e » G U TEE .(e*) and 
a e *-fi^e*)(* e .)GU^( e *)} 
={((x i ,p Q ;z Q ,u Q ),(a; l ,p Q ;r Q ,u a ))|(a; I ,z Q ) G lt(a; 1 ), 

/ = z Q and (a; 1 , r Q + u Q + C^z' 3 ) G IT^)}. 

Remark 3.2. One of the advantages of working in the Lie algebroids setting is that 
we can construct a local basis {e Q } of sections of E as follows. We take a local basis 
{e a } of sections of the vector bundle tu ■ U — » Q and complete it to a basis {e a , e^} 
of local sections of E. In this way, we have coordinates (x l ,y a ) = (x l ,y a , y A ) on E. 
In this set of coordinates, the equations which define the subbundle II are y A = 0. 
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So, we can consider (x l , y a ) as local coordinates on It. Moreover, if {e a , e A } is the 
dual basis of {e a , e^} of E* , then {e A } is a local basis of sections of 11°. So, from 
the definition of U-je e * , we deduce that {y a ,y a ,? A } is a local basis of sections of 
U<je e « ->• E* and, if {T,y A ,y a ,yA} is the dual basis of {^ a , y A , 7 a , 7 A }, then 
{y A } is a local basis of U% BE ,. Therefore, a local representation for the almost 
Dirac structure Du is 

^u{x\p a ) ={((x\ Pa ; z a ,u a ), (x*,p a ; r ai v a )) \ z A = v A = 0, v a = z\ 

and r a = -u a - Q 1 ab p 1 z b }. 

o 

We have used the canonical symplectic section il E on 7 E E* together with a 
vector subbundle U C E to define the almost Dirac structure Du- However there 
is a dual version of the above construction in which the almost Dirac structure is 
defined by a Poisson structure on 7 E E* together with a vector subbundle U C E. 

Let (E, [•, •], p) be a Lie algebroid and U be a vector subbundle of E. Consider the 
projection tt 2 : (7 E E*)* ->■ J5 defined as tt 2 = pr x ° (O^)" 1 , where pr x : 7 E E* ->■ £ 
is the projection on the first factor. If we consider local coordinates as before, using 
(I3.6[) . we have that 

ir 2 (x\p a ;r ai v a ) = (x\v a ). 
Now, we define the induced vector subbundle U 7 e e , of (7 E E*)* by 

U* 7EE , = (tt 2 )- 1 ^). 

Note that U^ EE , = Q e (U 7 e e ,), from the definition of U 7 e e , (see (I3.1[) ). Locally, 
U 7EE , is given by 

U 7EE ,{x\ Pa ) = {(x\p a ;r a ,v a ) | (x\v a ) e 

The annihilator of \i 7EE , is given, for each e* & E* , by 

(li;- B£; ») ( e *) ={X e , G 7 E E* I a e » (X e ») = 0, for all a e * eU* 7EE 4e*)} 

={Xt x t iPa \ = (x l ,p a ; z a , u a ) z a = and (x\u a ) e lt°(x 1 )}. 

On the other hand, we introduce the section II of the vector bundle A 2 (7 E E*) — > 
E* defined by 

n(a,p) = n E ((n E )- 1 (a),(n E )- 1 (p)), 

for a, (3 <E (7 E E*)*. H is the algebraic Poisson structure on the vector bundle 
7 E E* — !> E* associated with the symplectic section Q E . Denote by (tn : (7 E E*)* — > 
7 E E* the vector bundles morphism given by 

)J n (a) = -i a tt, for a G (7 E E*)*. 

Note that jjri = 

Then, using the above notation, the induced almost Dirac structure T>u on 7 E E* 
is given, for e* € E* , by 

D u (e*) ={{X e ,,a e .) £ 7 E E* x (7 E E*)*\a e . e l^- B£ ,.(e*) and 

-ttn(e*)(a e .) e (11^,)° (e*)}, 

whose local representation is 

£> u (y,p Q ) ={((a; l ,p Q ;2: Q ,u Q ),(a; l ,p Q ;r Q ,t; a )) | {x\v a ) e U(x 4 ), 
v Q = z" and (a; 2 , r Q + u Q + C^z' 3 ) e 

which coincides with (|3.7|) . 
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3.2. Implicit Lagrangian systems on a Lie algebroid. In this section, an 
implicit Lagrangian system on a Lie algebroid E is defined in the context of the 
induced almost Dirac structure Du on 7 E E* . As we shall see, the notion of implicit 
Lagrangian systems that is developed here can handle systems with degenerate 
Lagrangians as well as systems with nonholonomic constraints. Another description 
to address these systems was recently presented by Grabowska and Grabowski 
in [TB], where they use the notion of a Lie algebroid as a double vector bundle 
morphism. 

Let L : E — > M. be a Lagrangian function on the Lie algebroid (E, [•,•], p). 

First of all, we will recall the definition of the Legendre transformation in the 
context of Lie algebroids. Given a Lagrangian function L : E — > K, one can consider 
the Poincare-Cartan 1-section associated with L, 9^ <E T((7 E E)*), which is given 

6 L {e){Z e ) = (d 7EE L(e))(S e (Z e )) = p T {S e {Z e )){L), 

for e S E and Z e 6 7f E, S : 7 E E — > 7 E E being the vertical endomorphism defined 
in (12.51) . So, the Legendre transformation associated with L is defined as the 
smooth map ¥L : E — > E* defined by 

¥L(e)(e') = 6 L (e)(Z), 

for e,e' S E, where Z e 7 E E such that pr-^Z) = e', pr x : 7 E E E being the 
canonical projection over the first factor. For more details see |27j . 

The map FL is well-defined and its local expression in fiber coordinates on E 
and E* is 

FL(^) = (^|i). 

Now, we consider the isomorphism Ae ■ 7 E E* — > (7 E E)* between the vector 
bundles p^ : 7 E E* ->■ E and (r T )* : (7 E E)* E introduced in [37] and whose 
local expression is 

A E {x\p a ]z a ,u a ) = (x\z a ;u a + Qlpp-fZ 13 ,p a ). (3.8) 

Then, we define the map j E ■ (7 E E)* (7 E E*)* as j E = ^l E °A E l which is an 
isomorphism between the vector bundles (r r )* : (7 E E)* — > E and pr^ : (7 E E*)* — > 
E. From (|3.6I) and (13.81) . we deduce that the local expression of this isomorphism 
is 

j E (x l ,y a ;s a ,w a ) = {x l ,w a ;-s a ,y a ). (3.9) 

Now, define a differential operator D acting on the Lagrangian L : E — > R, which 
we shall call the Dirac differential of L by 

DL : E (7 E E*)*, DL = j E ° d^ '' E L, 

where d7 E L is the differential of L on the Lie algebroid 7 E E which is a section of 
(r T )* : (7 E E)* -» E. 

Using (|2.ip . (|2.4p and p.9p . we conclude that DL is represented in local coordi- 
nates by 

/ Fit Fit \ 

D ^ Q ) = (^;-^y Q )- (3-io) 

Now, we have all the ingredients to define an implicit Lagrangian system on a 
Lie algebroid. 

Definition 3.3. Let L : E — > K be a given Lagrangian function (possibly degen- 
erate) on a Lie algebroid (E, [•,•], p) and 11 C E be a given vector subbundle of 
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t : E — > Q. Denote by T>u the induced almost Dirac structure on the Lie alge- 
broid prolongation 7 E E* that is given by iS.ty) and DL : E — > (7 E E*)* the Dirac 
differential of L. Let P = ¥L(li) C E* be the image of V. under the Legendre 
transformation. 

An implicit Lagrangian system is a triple (L,U,X), where X is a section of 
the Lie algebroid prolongation t t : 7 E E* — > E* defined at the points of P, together 
with the condition 

{X,DL)€T> U . 

In other words, as P = ¥L(\L) C E* , X can be seen as a section of7 E E* -> E® Q E* 
defined at the points of XL ©q P and thus, we require that for each point e € II and 
with e* = FL(e) S P, we have 

(X( e> e%DL(e))€-Du(e*). 

Definition 3.4. A solution curve of an implicit Lagrangian system (L,VL,X) is 
a curve (x(t),y(t)) E ll(x(t)), t\ < t < t2, such that ¥L(x(t),y(t)) is an integral 
curve of the vector field p T (X) on E* , p T being the anchor map of the Lie algebroid 

Remark 3.5. One can consider the map i E ; E — »■ E ©q E* defined as the direct 
sum of the identity map on E, id : E — > E, and the Legendre transformation 
¥L : E — > E* . Denote by X the submanifold of E ©q E* defined as the image of U 
under is- Thus, % is locally given by 

BT 

% = {{x\y a ,p a ) S E x x El | {x\y a ) G U{x% p a = — }. 

oy a 

Another way to define the submanifold % is the following. Consider the map 
P(T E E)* '■ {T E E)* — > E ®q E* defined as the direct sum of the maps (r T )* : 
(T E E)* -> E and r T * ° (Ae)- 1 : {7 E E)* ->• E*. Recall that (t t )* : {7 E E)* -)• E is 
the projection of the dual vector bundle of the Lie algebroid prolongation of E over 
the fibration r, A E : 7 E E* ->■ (7 E E)* is the vector bundle isomorphism defined 
in (I3.8[) and r T : 7 E E* — > E is the projection of the Lie algebroid prolongation 
of E over r*. If we consider local coordinates introduced in Section the map 
P(je e) , is given by 

P(7ZE)- {x\y a ,s a , w a ) = (x\ y a , w a ). (3.11) 

Note that when E is the standard Lie algebroid, that is, E = TQ, then this map is 
the map p T * T Q ■ T*TQ -> TQ ®q T*Q defined in [S] (see Section 4.10 in g3]). 

Then, we can construct the map %e between E and E ©q E* by the composition 

ie = P(7B E y o A E o (O^-'oDL : E -> (7 E E*)* -> 7 E E* -> (7 E E)* -> E ®q E* , 

T>L being the Dirac differential of the Lagrangian function L and Q E being the flat 
map defined by the canonical symplectic section fl e (see (|2.8D ) . 

From (|3.6|) . p.8p . (I3.10|) and (|3.11[) , we have that the local expression of %e is 

Then, the submanifold C S ©q i?* can be defined as X — i E (U)- 
Then, a solution of an implicit Lagrangian system (L, XL, X) may be equivalently 
defined to be a curve (x(t),y(t),p(t)), where t\ < t < t 2 , whose image lies in the 
submanifold X C E ©q E* and such that (x(t),p(t)) is an integral curve of p T (X) 
and such that 



(X(x(t),y(t),p(t)),-DL(x(t),y(t))) € D u {x{t),p{t)). 
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Locally, using the preceding notation, (12.71) . (13.71) and (I3.10[) . we deduce that a 
solution curve (x l (t), y a (t),p a {t)) for an implicit Lagrangian system {L, 11, X) must 
satisfy the following equations 

dL 

(x\y a ) eU(x'), x l = p l a y a , Pa = — , 

Remark 3.6. If we consider the local coordinates on E introduced in Remark 13.21 
the implicit Lagrangian equations reduce to 

y A = 0, x^p l a y a , P a =g^ and Pa = -Q^—y* + p^—. 



3.3. Conservation of energy. Define the generalized energy El ■ E®q E* — s- 

R by 

E L (x,e,e*) = (e*,e) - L(x,e), 
where (x, e) ell and (x, e*) € P. 

Proposition 3.7. Let (x(t),y(t)), t\ < t < £2, oe an integral curve of a given 
implicit Lagrangian system (L,U,X) on a Lie algebroid E. Then, the function 
EL(x(t),y(t),p(t)) is constant in time, where p(t) = ^(x(t),y(t)). 

Proof. We give the proof using local coordinates. Then, from the definition of the 
generalized energy El, we have that 

dE L dL t i dL 

—r~ = V Pa + y Pa - TT- X - -= V ■ 

dt y F y F dx l dy° y 

dL 

As p a (t) — 7: — (x l (t), yP(t)), we deduce that 
dy a 

dE L a . dL .. 

V Pa - -7T-X 1 



dt y fe* 



Now, as (x l (t) , y a (t)) satisfies the implicit Lagrangian equations (|3.12[) . we know 
that 



/ dT dT \ 

[x\y<*)&U{x% x* = p* a y a and (x\p a + 6^— / - p l a —) g U°(^). 

ax 



Moreover, as = — C^ Q , the term y a £ 1 a p-Q-^y^ = 0. So, we conclude that 

□ 
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4. Hamilton- Jacobi theory for implicit Lagrangian systems 

Let (E, [■, •], p) be a Lie algebroid over a manifold Q with projection r : E — » Q 
and [L, U, X) be an implicit Lagrangian system on E. 

Theorem 4.1. Let 7 : Q — » i? 0q _B* fee a section of the canonical projection 
v : E <Sq E* — > Q smc/i i/iai 

7(Q) c DC, (4.1) 

and 

d E {pr E ,°j) lUxU = 0. (4.2) 

Denote by a £ T(i?) ifte section a — pr x °I°pr E , °7, where p^ : 7 E E* — > E is 
the projection on the first factor and pr B » : E ®q E* — > E* the projection over the 
second component. Then, the following conditions are equivalent: 

(i) For every curve c : / — > Q in Q such that 

c(t) = p(a)(c(t)), for all t, (4.3) 

the curve 7°c is a solution of the implicit Lagrangian system (L,U, X). 

(ii) 7 satisfies the Hamilton Jacobi equation for implicit Lagrangian 
systems: 

d E (E L °j) £11°. (4.4) 

Proof. We give the proof using local coordinates. We consider local coordinates 
(x l ) on an open subset V of Q and a local basis {e^} of sections of E defined on V, 
then we have the corresponding local coordinates {x l ,y a ) on E. Denote by p l a and 
C^g the structure functions of the Lie algebroid E with respect to (x l ) and {e a }. 

Suppose that 7(2;*) = (x z , ^"(x^), j a (x : ')). Then, the condition (|4.1j) means that 

BT 

(x\rW))eU(x l ) and j a (x i ) = —(x i ,^)), (4.5) 
and the condition (|4.2I) can be written locally as 

S^=(S^ + ^ e ^> v ' (4 - 6) 

for all v, w G 11 given locally by v = v^ep and w = w s e$. 

If c(£) = (c'(<)), it is easy to prove that equation (|4.3[) can be rewritten in local 
coordinates as 

^i)-fWi)R(cW). (4.7) 

Using the hypothesis (14.1[) (see its local expression (|4.5[) L we also have that 
equation (|4.4I) is locally written as 

for all i> = w a e Q G U. 

(i) =^> (ii) Assume that (i) holds. Therefore 

[(c ! (t), 7 a (cW))eU( c W), 

c\t)= 1 a {c{t)) P Mt)) 1 

(t) + e s a0 (c(tm s (c(t))^(c(t)) 

-pmt))^-{c{t)^{c{t)))y a {c{t)) G U°(c(t)). 
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Then, using ([3~5]l and P~!?j) . we have that at c(t), for w = iu Q e Q (c(t)) G U(c(i)), 

Then, we conclude that (ii) holds (see (J3ZEJ)). 

(ii) => (i) Suppose that (ii) holds, that is, condition (14.81) is satisfied. Let c : / — > 
Q a curve such that c(t) = p(a)(c(t)). Then, we have that 

c i (t)=r(c(tM(c(t)). 

Moreover, from (14.5j) . we also know that 

(c*(t),r(c(t)))eU(c(t)) and % ( c (t)) = —(c(t),r(c(t))), (4.10) 

Moreover, using f)4.5[) and (|4.6I) . we deduce that at c(i), for all w = w a e a (c(t)) G 
U(c(i)), 

= (07V Q + 7^„V + - Pi^y (4.1D 

Finally, from f|4. 10[) and (14. 1 1[) . we conclude that j°c is an integral curve of X. 

□ 



5. Examples 

5.1. The case U = E. This is perhaps the simplest case in which one has no 
constraints but the Lagrangian may be degenerate. In this case, the induced almost 
Dirac structure T>u — De is given by 

D E (e*) = {(X e ,,a e ,) e 7*E* x (7f,E*)* \ a e , = ^ E {e*){X e .)}. 

So, locally the equations defining the almost Dirac structure in this case are 

v a = z a and r a — —u a — G^pP-yZ 13 , 

where X e * = (x\p a ; z a , u a ) and a e * = (x\p a ;r a ,v a ). Then, a curve (x l (t) , y a (t)) 
in E is a solution of the implicit Lagrangian system if and only if 

Pa = or 1 x =PaV Pa = Pa &? ~ ■ 

This means that in this case, the condition of an implicit Lagrangian system is 
equivalent to the Euler-Lagrange equations for L (see Equations (2.40) in |27|). 

In the usual formulation of Lagrangian systems on a Lie algebroid, one must 
to restrict to admissible curves on the Lie algebroid E, that is, curves c(t) in E 
such that (c(t),c(t)) G 7^ t) E or, locally, if c(t) = (x 4 (t), y a (t)) then x l = p^y 01 . 
Notice that integral curves of an implicit Lagrangian system automatically satisfy 
this condition. 

In this case we can reformulate the Theorem 14. II as follows. 
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Corollary 5.1. Let 7 : Q — > E be a section of the Lie algebroid t : E — > Q such 
that 

d E (WL°j) = 0. 

Denote by a G F(E) the section a = pr x ° X °¥L°j, where pr 1 : T E E* — > E is the 
projection on the first factor. Then, the following conditions are equivalent: 

(i) For every curve c : / — >■ Q in Q such that 

c(t) = p{a){c{t)), for all t, 

the curve 7°c is a solution of the implicit Lagrangian system (L,E,X), 
or equivalently, 7 ° c is a solution of the Euler- Lagrange equations for the 
Lagrangian L. 

(ii) 7 satisfies 

d E {e L o 1 )^0, 

where Sl ■ E — > R is the energy function associated with L (see (2.39) in 
[27] ) which is defined as e L = p T (A)(L) - L, A E T(7 E E) is the Euler 
section and p T is the anchor map of the Lie algebroid t t : < J E E — > E (see 
Examvle \2~T\) . 

This result can be viewed as the Lagrangian version of the Theorem 3.16 in [27] . 

5.2. The case E = TQ. Let E be the standard Lie algebroid ttq : TQ — > Q. In 
this case, the sections of this vector bundle can be identified with vector fields on 
Q, the Lie bracket of sections is just the usual Lie bracket of vector fields and the 
anchor map is the identity map id : TQ TQ. A vector subbundle U of TQ is 
just a distribution Aq on Q. 

Moreover, in this case, the Lie algebroid (T E E* 1 [•, -] r , p T ) is the standard 
Lie algebroid (TT*Q, [-,-], id). So, the lift of the vector subbundle II = Aq to 
J E E* = TT*Q is just the distribution A T » Q on T*Q defined as 

At*q = (T^q)- 1 (A q ), 

ttq : T*Q — > Q being the canonical projection. Moreover, VLe = Qtq is the 
canonical symplectic 2-form on T*Q. Then, the induced almost Dirac structure 
X>u defined in Theorem 13. II is given, for each point z E T*Q, by 

T>A Q ={{v z ,a z ) E T Z T*Q x T* z T*Q\v z E A T , Q (z) and 

a z -U> TQ {z){v z )EA° T , Q {z)}. 

This almost Dirac structure coincides with the induced almost Dirac structure 
introduced by Yoshimura and Marsden in Thus, if we apply the results of 

Section T3. 2 1 to this particular case we recover the formulation of implicit Lagrangian 
systems develop in [33"] . 

Moreover, applying the Theorem l4.1l to this particular case one recover the result 
develop in |25j for standard implicit Lagrangian systems. 

Example 5.2. We are going to consider a simple example: the case of Euler- 
Poincare reduction. In this case, we consider the particular case when the manifold 
Q is a Lie group G and the distribution Aq is just TQ = TG (that is, the cases 
when 11 = E and E = TG). Let L : TG — > K be a left-invariant Lagrangian. Then, 
we have that (g(t),v(t)) E T g ^G, t% < t < £2, is a solution curve of the implicit 
Lagrangian system (L, TG, X) if and only ii g(t) is a solution of the Euler-Lagrangc 
equations for L on G and g(t) = v(t), for t\ <t < t^. 

On the other hand, let g be the Lie algebra associated with G which is a Lie 
algebroid over a point. As, L is a left-invariant function, we can consider the reduced 
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Lagrangian Z : g — > R, Z = Li fl . Taking U = g, we have that a curve £(i) e g is 
a solution of the implicit Lagrangian system (l,g,Y) if and only if it is a solution 
of the Euler-Poincare equations on g. Moreover, as well known, g(t) satisfies the 
Euler-Lagrange equations for L on G if and only if £(t) = g(t)~ g(t) satisfies the 
Euler-Poincare equations on g. Then, we conclude that (g(t),v(t)) is a solution 
curve of the implicit Lagrangian system (L, TG, X) if and only if = g(t)~ g(t) 
is a solution curve of the implicit Lagrangian system (l,g,Y) and g(t) — v(t). 

5.3. Nonholonomic mechanics on Lie algebroids. Let (E, [•, •], p) be a Lie 

algebroid. Nonholonomic constraints on the Lie algebroids setting are given by a 
vector subbundle li of E. In [10], the authors introduced the notion of a nonholo- 
nomically constrained Lagrangian system on a Lie algebroid £ as a pair (L,U), 
where L : E — > R is a Lagrangian function on E and U is the constraint subbundle, 
that is, it is a vector subbundle of E. 

If we consider local coordinates as in Remark I3~2l then a solution curve {x l (t), 
y a (t)), t\ < t < t2, on E for a nonholonomic system must satisfy the differential 
equations (see Equations (3.7) in [TU]) 

'x l = p\y a , 

dt\dy a ) + dyy abV Pa dx i ' 

y = o. 

So, a nonholonomic system is locally represented by an implicit Lagrangian sys- 
tem (L,VL,X) together with the condition 

rlT 

since (x l (t),p a (t)) = ¥L(x l (t),y a (t)), where FL is the Legendre transformation. 

Example 5.3. Consider the situation of Example 15.21 but with a non-trivial left- 
invariant distribution on G. Then, we have Q = G a Lie group, L : TG — > R a 
left-invariant Lagrangian and U = Ac, Ac ^ TG and Ac ^ {0}, a left- invariant 
distribution on G, that is, a standard nonholonomic LL system on G. As we 
have proved in general, (g(t),v(t),p(t)), t\ < t < t2, is a solution curve of the 
implicit Lagrangian system (L,TG, X) if and only if (g(t),v(t)) is a solution of 
the Lagrange-d'Alembert equations for L on G and p(t) = dL/dv(g(t),v(t)), for 

h<t< t 2 . 

On the other hand, this type of nonholonomic system on G may be reduced 
to a nonholonomic system on the Lie algebra g associated with G. The reduced 
Lagrangian Z : g — > R is Z = Li a and the vector subspace of g is given by 
c) = Ag(c). Then, one has a constrained system (Z,Z>) on g. So, a curve £(t) G g is 
a solution of the implicit Lagrangian system (Z, g, Y") if and only if it is a solution of 
the constrained Euler-Poincare equations (or the so-called Euler-Poincare-Suslov 
equations, see [13]) on g. 

6. Conclusion and Future work 

In this paper, we introduced the notion of an induced almost Dirac structure, 
and show how it leads to implicit Lagrangian systems on Lie algebroids. This 
provides a generalization of Lagrangian mechanics on Lie algebroids that can ad- 
dress degenerate Lagrangians as well as holonomic and nonholonomic constraints. 
Furthermore, we have obtained a Hamilton-Jacobi theory for such systems. 
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In future research, we aim to study the possibility of obtaining a Hamilton-Jacobi 
equation, as in the general case described in [55] , using the notion of Dirac algebroids 
given in |16j . In this case, the theory will include all important cases of Lagrangian 
and Hamiltonian systems, including systems with and without constraints, and 
autonomous and non-autonomous systems. 

Another interesting direction would be to generalize to Dirac algebroids the work 
in 25 that relates the Hamilton Jacobi theory of a Dirac mechanical system with 
symmetry and the Hamilton-Jacobi theory of the associated reduced Dirac system. 
Furthermore, the relationship between the various Hamilton-Jacobi theories for 
reduction of Dirac mechanical systems formulated on Dirac, Courant, and Lie alge- 
broids, and the formulations based on Lagrange-Poincare bundles |45) also remains 
to be studied. 



Acknowledgements 



This material is based upon work supported by the National Science Founda- 
tion under the applied mathematics grant DMS-0726263, the Faculty Early Ca- 
reer Development (CAREER) award DMS-1010687, the FRG award DMS-1065972, 
MICINN (Spain) grants MTM2009- 13383 and MTM2009-08166-E, and the projects 
of the Canary government SOLSUBC200801000238 and ProID20100210. We also 
want to thank the two anonymous referees for their helpful comments and sugges- 
tions. 



References 

R. Abraham, J. E. Marsden: Foundations of Mechanics, Addison- Wesley, 2nd edition, 1978. 
P. Balseiro, J.C. Marrero, D. Martin de Diego, E. Padron: A unified framework for mechanics. 
Hamilton-Jacobi equation and applications, Nonlinearity, 23 (2010), 1887—1918. 
L. Bates, J. Sniatycki: Nonholonomic reduction, Rep. Math. Phys., 32 (1) (1992), 99-115. 
A.M. Bloch : Nonholonomic Mechanics and Control, Interdisciplinary Applied Mathematics 
Series 24, Springer- Verlag New- York, 2003. 

A.M. Bloch, P.E. Crouch: Representations of Dirac structures on vector spaces and non-linear 
L-C circuits, in: Differential Geometry and Control (Boulder, CO, 1997), in: Proc. Sympos. 
Pure Math., vol. 64, AMS, Providence, RI (1997), 103-117. 

A. M. Bloch, P.S. Krishnaprasad, J. E. Marsden, R,. M. Murray: Nonholonomic Mechanical 
Systems with Symmetry, Arch. Rational Mech. Anal, 136 (1996), 21-99. 
J. F. Carinena, X. Gracia, G. Marmo, E. Martinez, M. Munoz Lecanda, N. Roman-Roy: 
Geometric Hamilton-Jacobi theory, International Journal of Geometric Methods in Modern 
Physics, 3(7) (2006), 1417-1458. 

J. F. Carinena, X. Gracia, G. Marmo, E. Martinez, M. C. Munoz Lecanda, N. Roman- 
Roy: Geometric Hamilton-Jacobi theory for nonholonomic dynamical systems, International 
Journal of Geometric Methods in Modern Physics, 7 (3) (2010), 431-454. 
J. Cortes, E. Martinez: Mechanical control systems on Lie algebroids, IMA J. Math. Control. 
Inform., 21 (2004), 457-492. 

J. Cortes, M. de Leon, J.C. Marrero, E. Martinez: Nonholonomic Lagrangian systems on Lie 
algebroids, Discrete and Continuous Dynamical Systems - Series A, 24 (2) (2009), 213-271. 
T.J. Courant: Dirac manifolds, Trans. Amcr. Math. Soc, 319 (1990), 631-661. 
T.J. Courant, A. Wcinstein: Beyond Poisson structures, Travaux en Cours, 27, Hermann 
Paris (1988), 39-49. 

M. Dalsmo, A.J. van der Schaft: On representations and integrability of mathematical struc- 
tures in energy-conserving physical systems, SIAM J. Control Optim., 37 (1) (1998), 54—91. 
Y.N. Fedorov, D. Zenkov: Discrete nonholonomic LL systems on Lie groups, Nonlinearity, 
18 (5) (2005), 2211-2241. 

M. Gotay, J. M. Nester, G. Hinds: Presymplectic manifolds and the Dirac-Bergmann theory 
of constraints, J. Math. Phys., 19 (1978), no. 11, 2388-2399. 

K. Grabowska, J. Grabowski: Dirac algebroids in Lagrangian and Hamiltonian mechanics, 
Journal of Geometry and Physics, 61 (2011), 2233-2253. 

K. Grabowska, P. Urbanski, J. Grabowski: Geometrical mechanics on algebroids, Int. J. 
Geom. Methods Mod. Phys., 3 (3) (2006), 559-575. 



21 



[18] M. Hcnncaux, C. Teitelboim: Quantization of gauge systems, Princeton University Press, 
Princeton, N.J., 1992. 

[19] P.J. Higgins PJ, K. Mackenzie: Algebraic constructions in the category of Lie algebroids, J. 

of Algebra, 129 (1990), 194-230. 
[20] D. Iglcsias-Ponte, M. de Leon, D. Martin de Diego: Towards a Hamilton-Jacobi theory for 

nonholonomic mechanical systems, Journal of Physics A: Mathematical and Theoretical, 41 

(2008), 015205. 

[21] D. Iglesias, J.C. Marrero, D. Martin de Diego, D. Sosa: Singular Lagrangian systems and 
variational constrained mechanics on Lie algebroids, Dynamical Systems, 23:3 (2008), 351- 
397. 

[22] J. Koiller: Reduction of some classical nonholonomic systems with symmetry, Arch. Rational 

Mech. Anal, 118 (1992), 113-148. 
[23] W.S. Koon, J.E. Marsden: The Hamiltonian and Lagrangian approaches to the dynamics of 

nonholonomic systems, Rep. Math. Phys., 40 (1997), 21-62. 
[24] W.S. Koon, J.E. Marsden: Poisson reduction of nonholonomic mechanical systems with sym- 
metry, Rep. Math. Phys., 42 (1/2) (1998), 101-134. 
[25] M. Leok, T. Ohsawa, D. Sosa: Hamilton-Jacobi Theory for Degenerate Lagrangian Systems 

with Holonomic and Nonholonomic Constraints, Journal of Mathematical Physics, 53 (7) 

(2012), 072905 (29 pages). 
[26] M. de Leon, J. C. Marrero, D. Martin de Diego: Linear almost Poisson structures and 

Hamilton-Jacobi equation. Applications to nonholonomic mechanics, Journal of Geometric 

Mechanics, 2 (2) (2010), 159-198. 
[27] M. de Leon, J.C. Marrero, E. Martinez: Lagrangian submanifolds and dynamics on Lie 

algebroids, J. Phys. A: Math. Gen., 38 (2005), R241-R308. 
[28] K.C.H. Mackenzie: General theory of Lie groupoids and Lie algebroids, London Mathematical 

Society Lecture Note Series, 213. Cambridge University Press, Cambridge, 2005. 
[29] J.E. Marsden, T.S. Ratiu: Introduction to Mechanics and Symmetry, second ed., in: Texts 

in Applied Mathematics, vol. 17, Springer- Verlag (1999). 
[30] E. Martinez: Lagrangian Mechanics on Lie Algebroids, Acta AppJ. Math., 67 (2001), 295-320. 
[31] E. Martinez: Geometric formulation of Mechanics on Lie algebroids, In Proceedings of the 

VIII Fall Workshop on Geometry and Physics, Medina del Campo, 1999, Publicacioncs de la 

RSME, 2 (2001), 209-222. 
[32] E. Martinez: Classical field theory on Lie algebroids: multisymplcctic formalism, preprint, 

arXiv:math.DG/0411352. 
[33] B.M. Maschke, A.J. van der Schaft, P.C. Breedveld: An intrinsic Hamiltonian formulation of 

the dynamics of L-C circuits, IEEE Trans. Circuits Syst., 42 (2) (1995), 73-82. 
[34] T. Mestdag, B. Langerock: A Lie algebroid framework for nonholonomic systems, J. Phys. 

A: Math. Gen., 38 (2005), 1097-1111. 
[35] J. Neimark, N. Fufaev: Dynamics on Nonholonomic systems, Translation of Mathematics 

Monographs, 33, AMS, Providence, RI (1972). 
[36] T. Ohsawa, A. M. Bloch: Nonholonomic Hamilton-Jacobi equation and intcgrability, Journal 

of Geometric Mechanics, 1 (4) (2009), 461-481. 
[37] T. Ohsawa, O. E. Fernandez, A. M. Bloch, D. V. Zenkov: Nonholonomic Hamilton-Jacobi 

theory via Chaplygin Hamiltonization, Journal of Geometry and Physics, 61 (8) (2011), 

1263-1291. 

[38] A.J. Van Der Schaft, B.M. Maschke: On the Hamiltonian formulation of non-holonomic 

mechanical systems, Rep. Math. Phys., 34 (1994), 225-233. 
[39] A.J. Van Der Schaft, B.M. Maschke: The Hamiltonian formulation of energy conserving 

physical systems with external ports, Arch. Elektr. Ubertrag., 49 (1995), 362-371. 
[40] A.M. Vershik, L.D. Fadeev: Lagragian mechanics in invariant form, Sel. Math. Sov., 1 (1981), 

339-350. 

[41] R.W. Weber: Hamiltonian systems with constraints and their meaning in mechanics, Arch. 
Ration. Mech. Anal, 91 (4) (1986), 309-335. 

[42] A. Weinstein: Lagrangian Mechanics and groupoids, Fields Inst. Comm., 7 (1996), 207-231. 

[43] H. Yoshimura, J.E. Marsden: Dirac structures in Lagrangian mechanics Part I: Implicit 
Lagrangian systems, Journal of Geometry and Physics, 57 (1) (2006), 133-156. 

[44] H. Yoshimura, J.E. Marsden: Dirac structures in Lagrangian mcchacnis. Part II: Variational 
structures, Journal of Geometry and Physics, 57 (2006), 209-250. 

[45] H. Yoshimura, J.E. Marsden: Dirac cotangent bundle reduction, Journal of Geometric Me- 
chanics, 1 (2009), 87-158. 



22 



M. LEOK AND D. SOSA 



M. Leok: Department of Mathematics, University of California, San Diego, 9500 
Gilman Drive, La Jolla, California, USA 

E-mail address: mleok@math.ucsd.edu 

D. Sosa: Departamento de Economi'a Aplicada y Unidad Asociada ULL-CSIC Geome- 
tri'a Diferencial y Mecanica Geometrica, Facultad de CC. EE. y Empresariales, Uni- 
versidad de La Laguna, La Laguna, Tenerife, Canary Islands, Spain 

E-mail address: dnsosaOull.es 



